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ABSTRACT 

Using a large sample of 26623 quasars with redshifts in the range 1.5 < z < 5.1 
with C IV A1549 A emission line in Fifth Data Release of the Sloan Digital Sky 
Survey (SDSS), we investigate the cosmological evolution of the Baldwin Effect, i.e. the 
relation between the equivalent width (EW) of the C iv emission line and continuum 
luminosity. We confirm the earlier result that there exists a strong correlation between 
the C IV EW and the continuum luminosity, and we find that, up to z 5, the slope of 
the Baldwin Effect seems to have no effect of cosmological evolution. A sub-sample of 
13960 quasars with broad C iv A1549 A emission line from SDSS is used to explore the 
origin of the Baldwin Effect. We find that C iv EW have a strong correlation with the 
mass of supermassive black hole (SMBH), and a weak correlation with the Eddington 
ratio, iBoi/iEdd- This suggests that the SMBH mass is probably the primary drive 
for the Baldwin Effect. 
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I INTRODUCTION 

The discovery of an anti-correlation between the equivalent 
width (EW) of C IV A 1549 emission line and the continuum 
luminosity measured at 1450A in quasar rest frame, known 
as the Baldwin Effect, was first made by Baldwin (1977) 
(e.g. see a review by Shields 2006). This correlation was 
explored in Active Galactic Nuclei (AGNs) for other broad 
emission lines, such as Lya, C ill A1909, He iv A1640, Mg 

II A 2800, Si IV A1400, O vi A1034, N v A1240, H/3 A4861; 
for narrow emission lines, such as [O ill] A5007, [O ll] A3727, 
[Nv] A3426, [Nlll] A3869; even for the Fe Ka emission line 
in X-ray spectrum (e.g. Wampler et al. 1984; Baldwin et al. 
1989; Netzer et al. 1992; Borson & Green 1992; Iwasawa & 
Taniguchi 1993; Laor et al. 1995; Francis et al. 1995; Green 
1996; Nandra et al. 1997; Groom et al. 2002; Dietrich et al. 
2002; Baskin & Laor 2004; Warner et al. 2004; Page et al. 
2004; Zhou et al. 2005; Netzer et al. 2006). Over the past 30 
years, a significant amount of effort has been expended to 
confirm this effect and explore its origin and evolution. 

The slope of line EW versus the luminosity for different 
ions depends on their ionization energy (e.g. Dietrich et al. 
2002). However we are not sure whether the slope of the 
Baldwin Effect evolves with redshift. 

Several interpretations about the origin of Baldwin ef- 
fect have been proposed. It is not clear whether the Baldwin 
effect is driven by the underlying physical parameters, such 
as the redshift. Eigenvector 1 of Boroson & Green (1992), 



the Eddington ratio (i.e. the ratio of the bolometric luminos- 
ity to the Eddington luminosity), the mass of supermassive 
black hole (SMBH), the gas metallically (e.g. Dietrich et 
al. 2002; Shang et al. 2003; Sulentic et al. 2007; Baskin & 
Laor, 2004; Bachev et al. 2004; Warner et al. 2004). Baskin & 
Laor (2004) found a strong correlation between the G iv EW 
and the Eddington ratio, and suggested that the Eddington 
ratio is the primary physical parameters driving the Bald- 
win effect. However, Shang et al. (2003) used the method 
of Spectral Principal Gomponent (SPG) to find there is no 
correlation between the Baldwin effect and the Eddington 
ratio or the mass as underlying physical parameters. 

With the reverberation mapping result, we may calcu- 
late the virial SMBH mass from the broad lines (e.g., H/3 
, Ha , Mg II , G IV ) (e.g. Kaspi et al. 2000; Bian & Zhao 
2004; Vestergaard & Peterson 2006; Salviander et al. 2007). 
It provides the possibility to explore these underlying phys- 
ical parameters for the origin of G iv Baldwin effect. 

In the past, most of the work was made by using rel- 
atively small samples. For a sample of 20 quasars, Baldwin 
(1977) first discovered this effect. Wills et al. (1999) used a 
complete sample of 22 QSOs from the Bright Quasar Survey 
(LBQS) to discuss the relation between the Baldwin effect 
and the Eigenvector 1. For a sample of 993 quasars from 
the Large Bright Quasar Survey (LBQS), Green, Forster 
& Kuraszkiewicz (2001) studied the relations of G iv EW 
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with the continuum luminosity and redshiftQ Dietrich et al. 
(2002) used a large sample of 744 type 1 AGNs {0 < z < 5) 
from International Ultraviolet Explorer (lUE) and Hubble 
Space Telescope (HST) to discuss the correlations of contin- 
uum and emission-line properties in the rest-frame ultravi- 
olet and optical spectral ranges. 

Here we use a large sample of 26623 quasars with C 
IV emission line from the Fifth Data Release (DR5, Adel- 
mann et al. 2007) catalog of the Sloan Digital Sky Survey 
(SDSS) to discuss C iv Baldwin Effect and its evolutionary 
effect, and we also use a large sample of 13960 quasars with 
broad C iv emission line to discuss the underlying physical 
parameters of the Baldwin Effect. The sample is described 
in §2, the results and the analysis are given in §3, and the 
conclusions are presented in §4. 



2 SAMPLE 

The sample used in this paper is selected from SDSS-I 
DR5 (Adelmann et al, 2007). SDSS-I was completed and 
DR5 (Adelman-McCarthy et al. 2007) was made public on 
30 June 2006. The SDSS Quasar Survey is continuing via 
the SDSS-II Legacy Survey. SDSS-I DR5 covers an imag- 
ing area of about 8000 square degrees and a spectroscopic 
area of about 5740 square degrees. SDSS DR5 catalog con- 
sists of 90,611 objects with automated spectral classification 
in SDSS DR5. Recently, Schneider et al. (2007) presented 
a quasars catalog from SDSS DR5, consisting of 77,429 
quasars with Mi < —22. The wavelength coverage of SDSS 
spectrum is 3800A to 9200A and we select 33943 quasars 
with redshift 1.5 < z < 6.0 to make C iv A 1549 available in 
SDSS spectrum. 

There are three stages during the pipeline of SDSS spec- 
tral fit. The first two stages are used to find lines to deter- 
mine the line redshift. The third stage is used to measure 
the line feature by a single Gaussian fit for the continuum 
subtracted spectrum. The chi-squared values can be used to 
evaluate the quality of the fit. □ We select the following pa- 
rameters of these quasars for our discussion: redshift (z), mjd 
of observation (mjd), plate ID (plate), fiber ID (fiber ID), a 
of fitted gaussian (sigma), error of a (sigma_Err), equivalent 
width (EW), error in equivalent width (EW_Err), continuum 
flux value at 1549A(/i549), of fit (chisq), and degree of 
freedom (nu, u). 

Through the SQL search language 0, we get 27545 ob- 
jects from SpecLine table in the SDSS DR5 with automated 
spectral classification of quasar, redshift 1.5 < z < 6.0, and 
the redshift confidence not less than 95%. 

922 objects with negative values of C iv EW or with 
EW_Err values larger than the values of its EW are excluded. 
At last, we have a sample of 26623 quasars with C iv emis- 
sion line to discuss the C iv Baldwin Effect and its evo- 
lutionary effect. In Figure 1 we show the histogram of the 
redshift. Only two quasars have redshift larger than five, 
their reshifts are 5.032. 

In order to explore whether the Eddington ratio or/and 

^ Their luminosity at 2500 Ais derived from the Bj photometric 
magnitude, since the LBQS spectra are not spectrophotometric 
http:/ /cas. sdss.org/astrodr5/en/help/docs/algorithm. asp 
http:/ /cas. sdss.org/astro/en/help/docs/sql_help. asp 
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Figure 1. Redshift distribution for the 26623 quasars. 

the SMBH mass is/are the underlying physical parameters 
for the origin of the Baldwin Effect, we calculate the SMBH 
mass for quasars with broad C iv line. For a sample of 130 
AGNs with HST archival spectra, Sulentic et al. (2007) sug- 
gested subtracting a narrow component of C iv A1549 A with 
FWHM < 1500km s"^ to discuss C iv A1549 A as an Eigen- 
vector 1 parameter. Then in our paper, from SDSS DR5, a 
sub-sample of 13960 quasars with broad C iv emission line 
is selected by the following criterions: (1)ct/(1 + z) > 3.3A , 
corresponding 1500km s"^ of C iv FWHM, (2)^7*^ < 1-5, 
showing acceptable one gaussian fit of C iv emission line. 



3 RESULTS AND DISCUSSION 

3.1 The slope and correlation coefficient of the 
Baldwin Effect 

The C IV EW is from the result of SQL searching. The con- 
tinuum luminosity La (1549 A) at rest frame is calculate from 
its flux (/1549) by the following formula, 

La(1549A) =47r4/i549 (1) 
where dL is the luminosity distance. For flat universe, 

di^^c{l + z)H^^ [ dz[{l + zf{l + Q,Mz)-z{2 + z)nK]~^'''} (2) 
Jo 

where c is the speed of light. Ho is the present Hubble con- 
stant, Q,M is the mass density, S7a is the vacuum density. The 
cosmological constants we use in our paper are — 0.7, 
= 0.3, Ho = 71km s"\ 

The following formula is used to express the relation 
between C iv EW and the continuum luminosity, 

log EW(C IV )= QH-/3 log La (I549A) (3) 

Figure 2 gives the correlation between C iv EW and 
continuum luminosity at 1549A for 26623 quasars from 
SDSS DR5. Taking the average of the upper and lower er- 
rors of C IV EW as the weight, we use the weighted least- 
square linear regression to parameterize the strong corre- 
lation (Press et al. 1992, P. 655). The correlation coeffi- 
cient is i? = —0.265 and the null hypothesis is less then 
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Figure 2. C IV EW (in units of A) versus tlie continuum luminos- 
ity at 1549 A(in erg s~i A~i) for total 26623 quasars. The solid 
line is the least-square linear fit eonsidering the error of C IV EW. 
The dash line is the least-square linear fit without considering the 
error of C IV EW. 



10"". The slope is /3 = -0.156 ± 0.0002 and the intercept 
is a = 8.549 ±0.011. An unweighted least-square fitting is 
also given (denoted with dash line in Figure 2). We notice 
that the centroid of the distribution of the data is between 
these two fit lines. The dash line is shifted vertically with 
respect to the solid line, which is due to that the C iv lines 
with lower EWs usually have the larger errors in logarithm 
space, and thus have smaller weights in our fitting. In the 
hereafter linear fittings, we just show the results of weighted 
least-square linear fittings, taking the mean value of the up- 
per and lower errors of C iv EWs as the fitting weights. 

By lUE spectra, Kinney et al. (1990) found that the 
slope is -0.17 ± 0.04 and Espey & Andreadis (1999) found 
it is —0.17 ±0.03. By the optical selected sample, Zamorani 
et al (1992) found it is -0.13 ± 0.03. For a sample of 158 
AGNs observed with the Faint Object Spectrograph on the 
Hubble Space Telescope, Kuraszkiewicz et al. (2002) found 
it is -0.198 ± 0.036. Dietrich et al. (2002) found the slope 
is -0.14 ± 0.02 for their total sample, and -0.20 ± 0.03 for 
objects with log ALa(1450A) > 44 erg . We get a slope, 
P = — 0.156±0.0002, which is consistent with some previous 
studies (Osmer et al. 1994; Laor et al. 1995; and Green 1996; 
Dietrich et al. 2002). In table 1, we list the results of the 
above researches. 

The redshift coverage of our sample is from 1.5 to 5.1 
(see Figure 1). We divide our sample into four parts accord- 
ing to the redshift, and the redshift bins are: 1.5 < z < 2, 
2 < z < 2.5, 2.5 < z < 3.5, and 3.5 < z < 5.1. Figure 3 
presents the C iv Baldwin Effect in different redshift bins. 
For each redshift bins, considering the error of C iv EW, we 
use the least-square linear regression to find the correlation 
between C iv EW and the continuum luminosity at 1549A. 
We find a strong correlation for different redshift bins. The 
correlation coefficients and the slopes are listed in Table 2. 



Steeper slopes are found in these four redshift-bins, respect 
to the total sample, and this is similar to the results of the 
paper from Dietrich et al. (2002). For these four bins, the 
slopes are almost the same, which shows no cosmological 
evolution of the C iv Baldwin Effect up to z ~ 5. 



3.2 The origin of Baldwin Effect 

The physical origin for the Baldwin Effect is still an open 
question. It is suggested that the Baldwin Effect is driven 
by the Eddington ratio and/or the SMBH mass (e.g. Baskin 
& Laor 2004; Warner et al. 2004). Using 81 BQS quasars, 
Baskin & Laor (2004) found that the correlation coefficient 
between C iv EW and the continuum luminosity is —0.154 
, and the correlation coefficient between C iv EW and the 
Eddington ratio is —0.581. They suggested that Baldwin 
Effect is driven by Eddington ratio, I/boI /^Edd ■ 

We calculate the SMBH mass from FWHM of C iv emis- 
sion line by the following formula (e.g. Vestergaard & Pe- 
terson 2006), 



log Mbh(C IV ) 



log 



FWHM(C IV ) 
1000 km s~i 



f ALa(1350 A) 
i lO-i^erg s-i 



0.53 



(6.66 ±0.01). 



where FWHM(C iv ) = 2.35o-/(l +z) in units of km s~^ , z 
is the redshift. The uncertainty of the SMBH mass is due to 
the errors of FWHM(C iv ), ALa(1350A), and the system 
error in equations (4) is from the uncertainties of the BLRs 
geometry and dynamics. The uncertainty of our calculated 
SMBH mass is about 0.5 dex (Vestergaard & Peterson 2006). 

We also calculate the Eddington ratio, i.g. the ratio of 
the bolometric luminosity (LboI ) to the Eddington lumi- 
nosity (I/Edd ), where LEdd = 1.26 x 1O^*(MBH/Af0)erg s"^ 
The bolometric luminosity is calculated from the monochro- 
matic luminosity at 1350A , LboI = c x L1350 , where L1350 
is the continuum luminosity at 1350A and c = 3 ~ 5 (Ric- 
ahrads et al. 2006), and we take c = 4. Considering the 
uncertainty of our calculated SMBH mass is about 0.5 dex, 
the uncertainty of the Eddington ratio is about 0.5 dex or 
more. 

In order to calculate the SMBH mass, we use a sub- 
sample of 13960 quasars with broad C iv emission line from 
SDSS DR5, which is stated in section 2. In Figure 4, we 
show the relations between the C iv EW with SMBH mass, 
the continuum luminosity at 1549 A , and Eddington ratio, 
Lboi/ LEdd- The results of these three relations are listed in 
Table 2. 

Our results show that the correlation coefficient for the 
relation between C iv EW and the continuum luminosity 
is larger than the other two relations. And the correlation 
coefficient of C iv EW with SMBH mass is larger than that 
of C IV EW with Eddington ratio. Therefore, by the larger 
sample, we can't confirm the suggestion that the Eddington 
ratio is the underlying physical parameter for C iv Bald- 
win Effect, and maybe the SMBH mass is a better driving 
parameter. 

Here we also present the relations between the C iv 
EW and the continuum luminosity, the SMBH mass, and 
ifioi/iEdd in different redshift bins, same to that in section 
3.1 (see Table 3). 



(4) 
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Figure 3. log-log correlation between the C IV EW and the continuum luminosity at 1549A (in erg s ^ A ^) for the DR5 sample of 
26623 quasars in different redshift bins. Each redshift bin is marked in every graph. Solid lines are the linear fits. 



Table 1. Results of the previous studies and of our work. Col.(l): the article; Col. (2): the sample of the work; Col. (3); redshift; Col. (4): 
the slopes of the linear fit for the Baldwin Effect of the studies. 



article 


sample 


redshift 


slope 


(1) 


(2) 


(3) 


(4) 


our paper 


26623 quasars from SDSS DR5 


1.5 < 2 < 5.1 


-0.156 ±0.0002 


Kinney et al.(1990) 


101 quasars and 88 Seyfert 1 galaxies 


z < 2.0 


-0.17 ±0.04 


Espey et al.(1999) 


200 AGNs 


z > 1.0 


-0.17 ±0.03 


Zamorani et al.(1992) 


five 'complete' samples 




-0.13 ± 0.03 


Kuraszkiewicz et al.(2002) 


158 AGNs 


0.001 < 2 < 3.822 


-0.198 ±0.036 


Dietrich et al.(2002) 


744 AGNs 


< z < 5.0 


-0.14 ±0.02 
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Table 2. Results for the samples in different redshift bins. Col.(l): different redshift bins; Col. (2): number of quasars in different redshift 
bins; Col. (3): correlation coefficient of the linear fit; Col. (4): slope of the linear fit together with its error; Col. (5): intercept of the linear 
fit with its error at logLx{15i9A) = 10^^ erg s~^ ; Col. (6): the possibility of null hypothesis. 



z 






N 


R 


slope 


intercept 


P 




(1) 






(2) 


(3) 


(4) 


(5) 


(6) 




1.5 < 


2 


< 2.0 


14370 


-0.360 


-0.231 ±0.0003 


1.980 ± 0.022 


< 10" 


4 


2.0 < 


Z 


< 2.5 


6076 


-0.349 


-0.251 ±0.0005 


2.050 ± 0.032 


< 10- 


4 


2.5 < 


Z 


< 3.5 


4338 


-0.384 


-0.251 ±0.0006 


2.176 ± 0.035 


< 10- 


-4 


3.5 < 


Z 


< 5.1 


1839 


-0.337 


-0.237 ±0.0015 


2.210 ±0.092 


< 10- 


4 



Table 3. Results for different relations for the 13960 quasars with broad C IV emission line. Col.(l): different correlations; Col. (2) 
correlation coefficient of the linear fit for the different correlations; Col. (3): slope of the linear fit together with the error; Col. (4) 
intercepts of the linear fit with the errors at logL;^(1549A) = lO^'^ erg s^^ , logMBH/^Q = 9, logLBoi/LEdd = ~li respectively; Col. (5) 
the possibility of null hypothesis. 



relations 






R 


slope 


intercept 


P 


(1) 






(2) 


(3) 


(4) 


(5) 


log C IV EW - 


log 


Lx (1549A) 


-0.372 


-0.187 ±0.0004 


1.954 ±0.023 


< 10 


log C IV EW - 


log 


Mbh/Mq 


-0.344 


-0.198 ±0.0005 


1.853 ±0.006 


< 10 


log C IV EW - 


log 


-^'BolZ-^'Edd 


-0.128 


-0.100 ±0.0006 


1.398 ±0.001 


< 10 -1 



4 CONCLUTIONS 

Baldwin Effect has been explored for a long time and we are 
still doing our efforts to explore it. In our paper, we discuss 
the Baldwin Effect and its evolutionary effect from a large 
sample of 26623 quasars (1.5 < z < 5.1) from DR5. And 
also, using a sample of 13960 quasars with broad C IV emis- 
sion line from DR5, we explore the physical parameters that 
drive the C iv Baldwin Effect. The main conclusions can be 
summarized as follows: (1) Baldwin Effect exists in the large 
sample of 26623 quasars from DR5 of SDSS; (2) According 
to the slopes in our four redshift-bins, up to 2 ~ 5, there's no 
evolutionary effect for the slope of Baldwin Effect, however, 
these slopes are relatively steeper than the slope for the to- 
tal sample of 26623 quasars; (3) By the 13960 quasars with 
SMBH masses estimate, we find that the relation between 
C IV EW and SMBH mass is stronger than that between C 
IV EW and Leoi/iEdd (from R = -0.344 to i? = -0.128), 
SMBH mass seems to be a better underlying physical pa- 
rameter. 
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Figure 4. the log-log correlation between the C IV EW with SMBH mass, the Eddington ratio, I/BoI /i'Edd i Eind the continuum luminosity 
for the 13960 DR5 sample. The solid line is the linear fit. 
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Table 4. Results of the samples in different redshift bins for the relations between the C IV EW with the continuum luminosity, the 
SMBH mass, and I/Bol/^Edd- Col.{l): different relations; Col. (2): different redshift bins; Col. (3): number of quasars in different redshift 
bins; Col. (4): correlation coefficient of the linear fit; Col. (5): slope of the linear fit together with its error; Col. (6): intercept of the linear 
fit with its error at logL_),(1549A) = W^^ erg s~^ , logMan/M© = 9, logLBoj/Lg^j = —1, respectively; Col. (7): the possibility of null 
hypothesis. 



relations 




z 


number 


R 


slope 


intercept 


P 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


log C IV EW - 


log Lx (1549A) 


1.5 < 2 < 2.0 


7981 


-0.484 


-0.244 ±0.0005 


1.936 ± 0.032 


< 10 -4 






2.0 < z < 2.5 


3227 


-0.452 


-0.234 ± 0.0008 


1.987 ±0.048 


< 10 -4 






2.5 < z < 3.5 


1855 


-0.513 


-0.251 ±0.0011 


2.103 ± 0.066 


< 10 -4 






3.5 < 2 < 5.1 


897 


-0.367 


-0.228 ±0.0025 


2.167 ±0.151 


< 10 -4 


log C IV EW - 


log Mbh/Mq 


1.5 < 2 < 2.0 


7891 


-0.354 


-0.207 ±0.0006 


1.823 ± 0.008 


< 10 -4 






2.0 < z < 2.5 


3227 


-0.412 


-0.234 ±0.0009 


1.855 ± 0.011 


< 10 -4 






2.5 < 2 < 3.5 


1855 


-0.413 


-0.212 ±0.0011 


1.943 ± 0.014 


< 10 -4 






3.5 < 2 < 5.1 


897 


-0.294 


-0.142 ±0.0020 


2.015 ± 0.025 


< 10 -4 


log C IV EW - 


log iBolZ-f'Edd 


1.5 < 2 < 2.0 


7981 


-0.310 


-0.247 ±0.0009 


1.898 ± 0.001 


< 10 -4 






2.0 < 2 < 2.5 


3227 


-0.148 


-0.123 ±0.0013 


1.890 ±0.002 


< 10 -4 






2.5 < z < 3.5 


1855 


-0.178 


-0.138 ±0.0017 


2.004 ± 0.002 


< 10 -4 






3.5 < 2 < 5.1 


897 


0.016 


0.012 ± 0.0030 


2.037 ± 0.004 


0.7182 



